UNCLASSIFIED 


_ AD  NUMBER _ 

ADB017412 

LIMITATION  CHANGES 
TO: 

Approved  for  public  release;  distribution  is 
unlimited. 


FROM: 

Distribution  authorized  to  U.S.  Gov't,  agencies 
only;  Test  and  Evaluation;  DEC  1976.  Other 
requests  shall  be  referred  to  Naval  Weapons 
Center ,  China  Lake,  CA. 


_ AUTHORITY 

NWC  ltr  17  Jun  1977 


THIS  PAGE  IS  UNCLASSIFIED 


THIS  REPORT  HAS  BEEN  DELIMITED 

and  cleared  for  PUBLIC  release 
UNDER  DOD  DIRECTIVE  5200,20  AND 
NO  RESTRICTIONS  ARE  IMPOSED  UPON 
ITS  USE  AND  DISCLOSURE, 

DISTRIBUTION  STATEMENT  A 

APPROVED  for  public  RELEASE; 

distribution  unlimited, 


NWC  TP  5815 


CJT 


Probability  Distribution  Associated 
With  m  on  n  Aerial  Duels 


by 

Frank  C.  Reed 

Systems  Engineering  and  Synthesis  Division 
Weapons  Department 

DECEMBER  1976 


Distribution  limited  to  U.S.  Government  Agencies  onlyi  test  end  evaluation;  9  September  1976.  Other  requests  for  this 
document  must  be  referred  to  the  Naval  Weapons  Center. 


Naval  Weapons  Center 


fffMWC 


CHINA  LAKE,  CALIFORNIA  93555 


Naval  Weapons  Center 

AN  ACTIVITY  OF  THE  NAVAL  MATERIAL  COMMAND 


R.  G.  Freeman,  III,  RAdm.,  USN  . Commander 

G.  L.  Hollingsworth . Technical  Director 


FOREWORD 

This  interim  technical  report  describes  the  mathematical  derivation  of  a  methodology  for 
determining  survival  probabilities,  exchange  ratios,  weapon  usage,  and  other  measures  of  effectiveness  in  air 
combat  duels.  Data  inputs  for  the  methodology  can  be  from  a  variety  of  sources  including  digital,  manned, 
and  “real  life”  combat  simulations.  The  methodology  is  capable  of  handling  m  on  n  type  engagements. 

This  work  was  conducted  under  several  tasks  including:  AirTask  A03P-03P2/008B/3F32-31 1-000, 
Anti-Air  Requirements  Study,  The  Joint  Short  Range  Air-to-Air  Missile  (JSRAAM)  Project,  and 
independent  research  funds. 

The  methodology  described  in  this  report  is  currently  being  programmed  and  will  be  available  to 
continue  evaluation  of  existing  air  combat  data  as  well  as  being  available  for  post  analysis  of  the  planned 
AIMVAL  flight  tests. 

This  report  was  reviewed  for  technical  accuracy  by  Dr.  E.  A.  Fay. 


Released  by 
W.  B.  PORTER,  Head 
Weapons  Department 
9  September  1976 


Under  authority  of 
G.  L.  HOLLINGSWORTH 

Technical  Director 


NWC  Technical  Publication  5815 


Published  by . Technical  Information  Department 

Manuscript .  5362/MS  76-1 1 1 

Collation  . Cover,  27  leaves 

First  printing . 125  unnumbered  copies 


L — w„. .. . 


...  k  .  •>  . 


in  jujus  iijiiiiaui.i. . 

MiMiaiiiaitaUiMiis 


UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  flWian  D K«  Enfttd) _ 

REPORT  DOCUMENTATION  PAGE  I  befoI^Se™“or» 

REPORT  NUMBER _  |2.  GOVT  ACCESSION  NO.  3.  RECIPIENT'S  CATALOG  NUMBER 

NWC-TP-5815  I _ ^ _ _  _ 

TITLE  (tnd  Subtitle)'  //T  lj*  TYPE  Of  REPORT  4  PERIOD  COVER 


kl  REPORT  NUMBER 


Probability  Distribution  Associated  jTith 
m  on  n  Aerial  Duels  t 


7.  authorc*; 


J5  TYPE  OF  REPORT  4  PERIOD  COVERED 

Test  and  evaluation 


6.  PERFORMING  ORG.  REPORT  NUMBER 


I.  CONTRACT  OR  GRANT  NUMBERfi) 


i 


Frank  C./ Reed 


[9  PERFORMING  ORGANIZATION  NAME  ANO  ADDRESS 


Naval  Weapons  Center 
China  Lake,  CA  93555 


Ml.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 


Naval  Weapons  Center 
China  Lake,  CA  93555 


10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  «  WORK  UNIT  NUMBERS 

AirTask  A03P-03P2/008B/ 
3F32-311-000 


'/Y  I12-  -WgBQRT.  OAT..E-  - - - 

//j  December  1976  | 


13.  ^ 


BER  OF  PAGES 


14  MONITORING  AGENCY  NAME  4  ADDRESS^//  dllltfnl  from  Controlling  Ollict)  I  IS.  SECURITY  CLASS.  fof  >»!•  npert) 


(w  -r 


UNCLASSIFIED 


ISa.  DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 


16.  DISTRIBUTION  STATEMENT  (a I  thlt  Rtport) 


Distribution  limited  to  U.S.  Government  agencies  only;  test  and  evaluation; 
9  September  1976.  Other  requests  for  this  document  must  be  referred  to 
the  Naval  Weapons  Center. 


is.  supplementary  notes 


G 


1 19.  xYvviOROS(Contlnu*onr9V9r9m*ld9ifrt9C999mrywuild9ntltybyblocknumb9r) 


Duels 

Aerial  Combat 
Survival  Probabilities 


Markov  Chains 
Weapons  Expenditures 


20.  ABSTRACT  (Continue  on  r«v«r««  old*  It  nactnary  md  Identity  by  block  number) 

See  back  of  form. 


DD  1  JAN *73  1473  EDITION  OF  1  NOV  6S  IS  OBSOLETE 


S/N  0102-014*  6601  I 


'03  Oi  ? 


_ UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Wian  Data  Enlari 


r 


UHLLASSiriEU 


_l.lijH|TY  CLASSIFICATION  OF  This  PAGEdFhan  Data  E nffd) 


(U)  Probability  Distribution  Associated  With  m  on  n  Aerial  Duels,  by 
Frank  C.  Reed.  China  Lake,  Calif.,  Naval  Weapons  Center,  December  1976. 
52  pp.  (NWC  TP  5815,  publication  UNCLASSIFIED.) 

(UfThis  report  presents  algorithms  for  determining  the  probability 
distributions  of  survivors  and  weapons  expended  in  m  on  n  simulated  aerial 
combat.  Both  shoot-look-shoot  and  shoot-switch-shoot  firing  doctrines  are 
considered.  For  shoot-look-shoot  firing  doctrines,  a  method  of  bounding 
numerical  results  due  to  the  ,<ghost'rt*'effect  in  simulated  combat  is  presented. 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAOEfWian  Data  Bnfrtd) 


NWC  TP  5815 


CONTENTS 


Section  1.  Introduction  .  3 

Section  2.  A  Special  Algorithm  for  1  on  1  Duels  .  5 

2.1  Survival  Probabilities  for  Two  Duelists  .  5 

2.2  Weapons  Usage  .  9 

Section  3.  Survival  Probabilities  for  m  on  n  Duels  .  13 

3.1  Pessimistic/Optimistic  Firing  Doctrine 

and  m/a  on  n/8  Duels  .  13 

3.2  A  Markov  Chain  Formulation  for  Handling  m/a  on  n/8  Duels  ..  14 

Section  4.  Weapons  Usage  for  m  on  n  Duels  .  22 

4.1  Updating  Formulae  . .  22 

4.2  Joint  Probability  Distributions  of  Survivors  and 

Missiles  Used  in  an  m/a  on  n/8  Duel  .  23 

Section  5.  Shoot-Switch  Firing  Doctrines  .  25 

5.1  Bounds  on  the  Number  of  Non-zero  Probability  States .  25 

5.2  Shoot-Switch  With  a  =  0,  8  =  0  .  27 

5.3  Modified  Shoot-Switch,  a  -  0,  8  =  0  .  28 

5.4  Shoot-Switch,  a  >  0,  8  £  0  .  29 

Appendixes : 

A.  A  Numerical  Example  .  31 


B.  The  Numerical  Example  Continued  To  Include  Weapons  Usage  ...  43 


1 


NWC  TP  5815 


« 


Section  1.  INTRODUCTION 


This  report  is  written  to  derive  survival  probabilities,  exchange 
ratios,  weapons  usage,  and  other  measures  of  effectiveness  in  aerial 
duels  based  on  simulated  aerial  combat.  It  is  convenient  to  think  of  an 
m  on  n  duel  as  taking  place  between  two  sides.  Side  I  has  m  members  and 
side  II  has  n  members.  Simulated  aerial  combat  geometry  is  generated  by 
allowing  all  (m  +  n)  participants  to  fly  for  the  entire  duration  of  the 
duel. 


From  the  trajectories  on  these  runs,  each  combatant  is  assigned 
firing  opportunities  consistent  with  missile  flight  characteristics, 
envelope  considerations,  an  appropriate  firing  doctrine,  and  weapons 
loadout.  The  member  of  the  side  being  fired  at  is  assigned  based  on 
both  geometry  and  previous  firing  history.  Thus,  an  attempt  is  made  to 
minimize  overkill  and  spread  kill  potential  uniformly  over  the  history 
of  the  duel.  The  kill  probability  associated  with  the  launch  is  based 
on  relative  geometries,  single-shot  kill  probabilities  and  salvo  size. 

Unfortunately,  an  unavoidable  difficulty  is  induced  by  the  geometry. 
Whenever  m  >  1  or  n  >  1,  it  is  possible  for  a  combatant  who  would  have 
died  early  in  the  engagement  to  continue  to  draw  attention  and  fire  from 
an  adversary  who  would  otherwise  be  free  to  engage  viable  targets.  Thus, 
a  pursuer  has  the  chance  of  firing  at  "ghosts"  or  decoys,  rather  chan 
directing  his  attention  at  real  threats. 

In  addition  to  the  geometry  problem  there  is  apparently  some  dif¬ 
ficulty  in  making  probability  calculations  correctly.  Previous  attempts 
to  make  the  probability  calculations  using  a  recursive  algorithm  violate 
a  critical  independence  assumption.  A  clear  discussion  of  this  fallacy 
in  duel  calculations  is  contained  in  a  report  by  Lincicum.1  It  is  a 
common  error  made  by  a  number  of  contractors  and  government  facilities. 
This  report  considers  both  aspects  of  these  problems  for  general  m  on  n 
duels.  Section  2  provides  a  special  algorithm  which  correctly  obtains 
survival  probabilities  for  the  case  m  =  n  =  1.  Weapons  usage  is  also 
determined  for  this  case. 


*  Naval  Weapons  Center.  Determinations  of  Kill  Probabilities  and  Exchange  Ratios  for  Multiple  Firings  h\  Two 
Combatants,  by  L.  Lincicum.  China  Lake,  Calif.,  NWC,  April  1973.  (TN  3007-129,  publication  UNCLASSIFIED.) 
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Section  3  presents  a  different  algorithm  for  the  general  m  on  n  duel. 

The  concept  of  m/a  on  n/3  duels  is  introduced  as  a  method  of  handling  the 

ghost  problem  in  shoot- look-shoot  aerial  duels.  The  problem  is  then 
formulated  as  a  Markov  process  with  appropriate  updating  formulae  at 
state  space  transition  times.  Survival  probabilities  and  exchange  ratio 
formulae  are  also  presented. 

Section  4  extends  the  approach  discussed  in  Section  3  to  include 

weapons  usage.  This  extension  is  made  at  the  cost  of  what  appears  to 

be  an  exceedingly  large  state  space.  In  application,  the  actual  size  of 
the  state  space  one  must  consider  may  be  much  smaller. 

Section  5  discusses  shoot-switch  firing  policies  with  updating 
formulae  when  the  "ghost"  problem  is  ignored.  The  problem  of  combining 
m/a  on  n/8  duels  with  a  shoot-switch  firing  doctrine  is  discussed. 

Appendixes  A  and  B  illustrate  the  updating  formulae  in  Sections  3 
and  4,  respectively.  The  important  observation  in  this  example  is  that, 
even  though  the  potential  state  space  is  quite  large,  only  a  small 
fraction  of  these  states  are  ever  used.  Thus,  a  computer  program  which 
only  introduces  states  as  required  may  be  quite  efficient. 
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Section  2.  A  SPECIAL  ALGORITHM  FOR  1  ON  1  DUELS 


When  m  =  1  and  n  =  1,  we  are  concerned  with  an  important  special 
case  of  m  on  n  duels.  The  algorithm  presented  here  is  more  efficient 
than  the  general  algorithm  of  the  next  section  and  it  is  the  preferred 
method  of  analysis  of  1  on  1  duels.  The  "ghost"  problem  is  not  a  con¬ 
sideration  in  1  on  1  duels. 


2.1  SURVIVAL  PROBABILITIES  FOR  TWO  DUELISTS 


In  this  section,  which  is  based  on  a  report  by  Reed,2  we  describe 
an  algorithm  procedure  for  generating  time- dependent  survival  functions 
for  two  duelists  whom  we  call  I  and  II.  We  allow  time  delay  between 
launch  and  impact,  so  that  the  probability  both  duelists  are  killed  may 
be  greater  than  zero. 


We  assume  from  simulation  records  that  I  launches  his  jth  missile 
at  II  at  time  tj(I),  j  =  l,2,--*,r,  where  r  <,  Lj_,  the  weapon  load  for  I. 
We  assume  that  the  corresponding  intercept  times  Tj(I)  are  such  that 
T j ( I)  <  Tj+i(I)  for  all  j  and  we  allow  Pj(I),  the  kill  probability  of  I's 
jth  shot,  to  vary  with  j.  Similarly,  II  launches  his  ith  missile  at  I 
at  times  ti(II) ,  i  =  1,2, •••,s  £  L2,  II 's  weapon  loadout.  The  ith  impact 
of  II's  missiles  occurs  at  T^(II)  with  Ti(II)  <  Ti+i(II)  for  all  i  with 
probability  of  kill  Pi(II). 


We  define  P^j  as  the  probability  that  I  is  killed  at  time  1^(11)  and 
II  is  killed  at  time  Tj(I),  for  j  =  l,2,***,r;  i  =  l,2,**-,s.  We  define 


’Oj  as  the  probability  II  is  killed  at  Tr(I)  and  I  survives  the  duel 
Similarly,  P-^q  is  the  probability  I  is  killed  at  1^(11)  and  II  survives 
the  duel.  We  let  Pqq  be  the  probability  that  both  I  and  II  survive  the 
duel. 


The  key  to  constructing  a  valid  algorithm  comes  from  considering  a 
general  diagram  which  describes  the  interactions  between  missile  launches 
and  subsequent  impacts  of  the  two  duelists.  As  an  example,  consider  the 
following  figure: 


^  Naval  Weapons  Center.  Survival  Probabilities  for  Two  Duelists,  by  F.  C.  Reed.  China  Lake,  Calif.,  NWC,  June 
1972.  (IDP  3381,  publication  UNCLASSIFIED.) 
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We  say  that  the  duel  "ends"  when  one  or  the  other  duelist  is  killed. 
This  does  not  mean  that  the  first  duelist  killed  is  the  only  one  kille  . 
With  this  in  mind,  the  first  chance  for  the  duel  to  end  occurs  with 
the  first  impact.  In  the  figure  above  Tx(I)  is  the  earliest  impact  time. 


The  duel  "ends"  at  Ti(I)  with  probability  Pi(I).  Given  that  the 
duel  ends  with  II 's  death  on  this  first  impact,  it  is  a  simple  proba- 
bilistic  calculation  to  obtain  Pu  for  i  =  1,2,  ---.s.  One  merely  loo  s 
at  the  figure  to  see  which  launches  II  gets  off  before  he  is  killed. 

In  this  example,  II  gets  two  shots  off,  so  the  appropriate  conditional 

probabilities  are 


3  <  i  <  s 


P11|ll  killed  at  T^I)  -  PjUD 

P21 1 II  kllled  at  Tl^  “  ^  "  Pl*11^ 

P  | II  killed  at  TjU)  =  0 

P  |  II  killed  at  1^(1)  -  (1  -  P-^lDMl  -  P2(H)) 

The  unconditional  probabilities  P^i  for  0  <  i  ^  s  are  given  by 

P11  =  P1(I)  P1(II) 

P  =  P1(I)(1  -  P^ID)  P2(II) 


■pwintpm  w  ^ 
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Pn  =  0  i  =  3,4,  •  •  •  ,s 

p0i  =  p1(i)(i  -  P1(ii))d  -  p2(ii)) 


With  probability  (1  -  Pi(I))  the  duel  does  not  end  at  T^(I).  Thus 
with  probability  (1  -  P^(I))  one  gets  a  figure  like  that  above  with  the 

vector  ti(I)l^(I)  eliminated.  Proceeding  in  this  way  one  generates  all 
probabilities  required,  stopping  when  all  vectors  have  been  eliminated. 

We  now  turn  to  an  algorithm  which  formalizes  these  ideas.  We  let 

w1’  ■  vi(i)  -  wn>  -  wn)  ■ 1 


where 

T  >  max  [ T r ( I ) ,  Tg ( II)  ] 

We  also  let 

P0(I)  =  P0(H)  =  0 

For  any  ig,jg  we  proceed  as  follows: 

1.  If  ig  =  s+1  and  jg  =  r+1,  set 

s  r 

p0o  ■  n  <i  -  h(n»  n  o  -  pj<t» 

i-0  j-0 

and  end  the  algorithm;  otherwise,  go  to  2. 

2.  If  Tio(II)  <  x Jq (I) ,  go  to  3;  otherwise,  go  to  6. 

3.  Form 

S(I)  =  |j  :  tj  (I)  <  tj(I)  <  Tio<II>| 


If  S(I)  =  <p,  go  to  4;  otherwise,  go  to  5. 


and  go  to  1  with  ig  replaced  by  ig+1. 


Set 

*0  1  max  S(I) 

pi  o =  pi (n)  n  (i-mii»  n  a-*AD) 

0  0  i-0  j-0 

V1  j-i 

pi0j  =  pi  77  a  -  vid)  77  u  -  pk(D) 

i=0  k=0 

for  je  S ( I)  ,  and  go  to  step  1  with  ig  replaced  by  ig+1 
Form 

S(II)  =  I  i  :  t,  (II)  <  t.(II)  <  T.  (1)1 

I  x0  J0  I 

If  S(II)  =  4),  go  to  7;  otherwise,  go  to  8. 

Set 

V1  jo'1 

poi  =  pi  <*>  77  d  -  mid)  77  (i  -  md> 

J°  j°  i=o  3=0  j 


and  go  to  step  1  with  jg  replaced  by  jg+1. 
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8.  Set 


P  =  P 

OJn  : 


'0 


max  S ( II) 

a)  77  a 

i=0 


Jo"1 

p.dD)  n  (i 

j-0 


P.(D) 


iJn  "  J 


i-i  Jo"1 

p4  d)Pi(ii)  fj  (i  -  pkUD>  77  d  - 

k=0  j=0 


Pj(D) 


for  ie  S(II)  and  go  to  step  1  with  jg  replaced  by  jg+1. 
The  algorithm  begins  with  ig  =  jg  =  1. 


Survival  probabilities  for  the  two  duelists  are  easily  obtained  when 
the  P.jj  have  been  obtained  from  the  algorithm  above.  If  Sj(t)  is  the 

probability  I  survives  through  time  t,  we  have 

i  r 

w11”  - 1  -£  £  pkj 

k=l  j  =0 

and  I's  duel  survival  is 

s  r  r 

si  - 1  -  £  E  pkj  ■  E  poj 

k=l  j=0  j=0 

Similarly,  if  SII(t)  is  II's  probability  of  survival  through  t, 
we  have 


sII(Tj(D) 


l  - 


EE 

i=0  k=l 


ik 


9 


*tj"*i**r.  •*' 


n^yyjiYWtfj>.  v;  «fa  -  • .  *4 


NWC  TP  5815 


and  II 's  duel  survival  is 


SII  "  1  ‘  £  £ 


2  PiO 


i=0  k=l 


It  is  important  to  notice  that  the  above  algorithm  does  not  depend 
on  a  shoot-look-shoot  firing  doctrine.  It  only  depends  on  the  kth 
missile  fired  by  I(II)  not  reaching  I I ( I)  before  the  k+rth  missile  fired 
by  I (II)  reaches  II (I).  We  now  turn  to  weapons  usage. 


2.2  WEAPONS  USAGE 

We  define  as  the  random  variable  which  takes  on  value  1  if  I  is 
alive  and  0  if  I  is  dead.  The  random  variable  X2  is  defined  similarly 
for  II.  The  random  variable  takes  on  the  integer  values  associated 
with  I's  use  of  weapons.  The  random  variable  M2  is  defined  similarly 
for  II.  The  joint  distribution 


P[X.  =  x.  ,  X,,  =  x0 ,  M.  =  m.  ,  M„  =  m„]  =  P 

ll  1’  2  2’  1  1*  2  2  xlx2mlm2 


P  =  P 
llrs  00 


P  =  p 

10jL2  Oj 


j  =  1,2, • • • ,r 


P  =  p 

OlLji  iO 


i  =  1,2,  •  •  • , s 


pnm.  t.  1  pno  “  ^  prn  ^  p 


ool1l2 


00  Oj  iO 

j=l  J  i=l 


From  Px-^x2m^m2  we  may  compute  the  following  means,  variances,  and 


covariances : 


E(xi>  ■  £  po. 
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s 


E<V  ■  £  P10 

i=0 


Var(XL) 


Var(X2) 


s 


£  PiO 
i=0 


Cov(X1,X2) 


r 


E<M1>  ■  rP00  +  £  JP0j 

j=l 


Var^) 


e2(m1) 


Var(M2) 


i=l 


-  e2(m2) 
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Section  3.  SURVIVAL  PROBABILITIES  FOR  m  ON  n 
SHOOT-LOOK-SHOOT  DUELS 


We  now  want  to  consider  more  general  m  on  n  shoot -look-shoot  duels. 
For  this  report  shoot-look-shoot  implies  each  combatant  may  have  at  most 
one  missile  in  flight  at  a  time.  The  procedure  will  be  to  define  a  state 
space  for  which  the  Markov  assumption  holds  and  then  provide  the  proba¬ 
bilistic  updating  required.  First,  we  introduce  the  concept  of  m/a  on 
n/6  duels  as  a  means  for  handling  the  "ghost"  problem. 

3.1  PESSIMISTIC/OPTIMISTIC  FIRING  DOCTRINE 
AND  m /a  ON  n/0  DUELS 

From  the  simulated  aerial  engagement,  missile  simulation,  and  rules 
of  assignment,  one  obtains  a  time  history  t]_  <  t£  •  •  •  £  tjj  where  each  t^ 
is  either  a  missile  launch  or  impact.  For  each  tfc  the  combatant  firing 
the  missile  salvo,  the  intended  victim,  and  missile  salvo  kill  probability 
are  all  known.  The  victim  is  specified  by  the  geometry  of  the  engagement 
and  does  not  account  for  the  fact  that  he  may  not  be  alive.  This  means 
that  whenever  a  combatant  is  being  shot  at  for  the  ith  time  where  i  >  2 
and  kill  probabilities  are  greater  than  zero,  there  is  a  nonzero  proba¬ 
bility  that  the  shot  is  wasted.  We  think  of  such  a  constraint  as  a 
"pessimistic  firing  doctrine,"  since  it  is  pessimistic  from  the  point 
of  view  of  the  side  doing  the  firing. 

This  geometric  inaccuracy  affects  both  sides  and  one  would  like  to 
think  the  results  have  a  balancing  effect.  Intuitively  one  would  ex¬ 
pect  the  effect  to  be  more  serious  for  the  better  side.  Every  time  the 
better  side  kills  an  opponent  it  is  replaced  by  a  nonlethal  decoy.  These 
nonlethal  decoys  tend  to  increase  the  survivability  of  opponents  of  the 
better  side.  For  1  on  1  duels  the  "ghost"  effect  is  not  serious,  for  a 
nonlethal  decoy  does  not  degrade  the  survivability  of  the  winner. 

To  compensate  for  the  geometry  inadequacy,  we  consider  different 
target  allocations  than  those  implied  by  the  geometry.  Targets  are  as¬ 
signed  using  an  assumption  that  is  easy  to  implement  with  the  existing 
data  base.  One  merely  assumes  that  if,  in  the  original  geometry,  the 
target  is  dead  when  a  missile  is  launched  at  him,  then  this  firing  may  be 
converted  to  a  live  target,  if  there  is  one.  The  ability  to  always  con¬ 
vert  is  called  an  "optimistic  firing  doctrine"  and  intuitively  it  gives 
the  pursuer  more  capability  than  he  would  have  in  a  real  duel. 

We  expand  on  this  concept  by  defining  an  m/a  on  n/8  duel  as  an  m 
on  n  duel  in  which  side  I  has  optimistic  firings  with  probability  a, 
side  II  has  optimistic  firings  with  probability  8>  and  pessimistic  firings 
occur  with  the  complements  of  these  probabilities.  If  E(m/a,  n/8)  is  the 
measure  of  effectiveness  for  side  I  associated  with  an  m/a  on  n/8  duel 
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and  E(m,  n)  is  the  unknown  measure  of  effectiveness  associated  with  an 
m  on  n  duel  with  the  geometry  deficiency  removed,  then  intuitively, 

E(ra/0,  n/1)  <  E(m,  n)  <  E(m/1,  n/0) 

This  should  allow  one  to  bracket  the  payoff  to  side  I  with  two  duel  cal¬ 
culations.  Both  E(m/0,  n/0)  and  E(m/1,  n./l)  will  l^e  in  this  interval 
and  one  wonders  if  there  exists  an  ag  and  $q  so  that  E(m/otQ,  n/8o)  best 
estimates  E(m,  n)  over  a  wide  class  of  m  on  n  duels. 

3.2  A  MARKOV  CHAIN  FORMULATION  FOR  HANDLING 
m/a  ON  n/0  DUELS 

We  recall  that  the  basic  inputs  to  the  problem  are  times 
tl  ^  tz  <  •  •  •  <  tN  where  each  t^  is  either  an  impact  or  launch  time. 
Moreover,  for  each  t^  the  launcher,  target,  and  kill  probability  are 
known.  Our  goal  in  this  section  is  to  define  an  appropriate  state  space 
so  that  the  stochastic  process  defined  on  this  state  space  at  times  t^ 
satisfies  the  Markov  assumption.  Thus  one  may  update  probability  dis¬ 
tributions  at  time  tfc  by  knowing  the  distributions  at  time  t^-i  and  the 
appropriate  transition  probabilities. 

To  formulate  the  problem  as  a  Markov  chain,  one  must  be  careful  to 
construct  a  state  space  for  which  the  Markov  assumption  holds.  Once  this 
is  done  the  problem  is  merely  one  of  properly  updating  the  probability 
distributions  of  the  state  of  the  process  at  launch  and  impact  times. 

We  define 


X1  (xll’  X12 ’  '  "  ’  Xlm) 

where  xn  =  1  if  the  ith  member  of  side  I  is  alive  and  xn  =  0  if  the 
ith  member  of  side  I  is  dead.  The  components  of  vector 


X2  (x21’  X2n) 

are  defined  similarly  for  side  II.  We  also  have 
T1  =  (tll’  t12’  '  “  *  tlm) 

where  tj_i  =  0  if  the  ith  member  of  side  I  has  no  missile  in  flight  and 
tn  =  j  if  the  missile  of  the  ith  member  of  side  I  is  directed  at  the 
jth  member  of  side  II.  In  a  similar  way. 


T2  ~  (t21’ 


•’  t2n) 


where  t2j  =  0,  1,  •••,  m.  The  state  space  is  determined  by  the  vector 
(Xj_  9  X2 ,  Ti,  T2)  and  the  associated  Markov  process  has  2nr+-n(n  +  l)m(m  +  l)n 
states.  Clearly,  if  m  and  n  are  large,  the  number  of  states  is  large. 
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In  order  to  reduce  the  size  of  the  state  space  there  Is  a  real 
temptation  to  let  (X]_,  X2)  describe  the  state  space.  If,  however,  the 
missile  of  element  i  of  side  I  is  to  impact  element  j  of  side  II  at  a 
transition  time  and  if  xii  =  0,  we  do  not  know  whether  to  allow  the 
impact  or  not.  Element  i  of  side  I  may  or  may  not  have  been  alive  at 
the  time  of  the  missile  launch.  The  state  space  description  described 
here  is  slightly  more  generous  than  actually  required;  however,  this 
state  space  description  allows  us  to  assess  the  effect  of  the  geometry 
deficiency  described  previously. 

We  are  now  ready  to  consider  the  updating  of  the  probability  dis¬ 
tribution  defined  on  the  state  space  of  the  Markov  process  with  the 
possibility  of  optimistic  firings.  It  must  be  stressed  that,  with 
optimistic  firings,  targets  other  than  those  implied  by  the  geometry 
may  come  under  consideration.  The  state  space  we  have  constructed 
always  allows  us  to  decide  who  the  target  is,  even  if  optimistic  firings 
are  considered.  *- 


We  shall  make  use  of  the  following  notation: 

An  r-vector  of  l's 
0r  An  r-vector  of  0's 

6  An  r-vector  of  0's  with  the  exception  of  component 

ri  i  which  is  a  1 

6x  A  scalar  which  is  0  if  x  ■  0  and  1  if  x  >  0 

Pli(tk)  Kill  probability  of  salvo  launched  by  member  i 

of  side  I  and  impacting  at  t^ 

P„,(tk)  Kill  probability  of  salvo  launched  by  member  j 

J  of  side  II  and  impacting  at  t^ 

For  s  =  0,l,***,r-l,  we  define 


kr(i,  s)  =  i  +  s 

=  i  +  s  -  r 

With  respect  to  side  I  we  define  s^* 
for  which 


if  i  +  s  <  r 
if  i  +  s  >  r 

to  be  the  smallest  integer  s 


si*>  *  1 
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If 


then 


kma,  sL*) 


0 


With  respect  to  side  II,  we  define  s2*  to  be  the  smallest  integer  s 
for  which 


x2kn(j,  s2*) 


If 


then 


0 

n 


kn(j ,  s2*)  i  0 

We  define  P„  „  „  (t)  to  be  the  probability  that  the  stochastic 

x1x2t1t2 

process  under  consideration  is  in  state  Xi,X2,T^,T2  at  time  t.  Now 
for  all  tk  we  want  to  take  Px1x2T1T2(tk-l)  ’  the  Probability  distribution 
defined  on  the  state  space  at  tk_i  and  update  it  based  on  the  Markov 
property  to  obtain  Px2X2TiT2^tk^ '  assume  tbat  at  b0  <  bl 

PL  L  0  0  (t(P  =  1 

m  n  m  n 

Now  assume  that  at  tk  member  i  of  side  I  shoots  (according  to  the  geom¬ 
etry)  at  member  j  of  side  II.  Assuming  Px  X  T  T  ^k-l^  *s  ^nown  tben 
the  probability  updating  can  be  given  as  12  12 


f,WlT2(t|‘)  *  (X  -  SifWl 


T2^Ck-l^ 


"*■  Px  X  T  -16  T  ^bk-l^ 

xli  X2j  X1X2T1  j°miT2 


"*■  ^  ~  a^PX  X  T  -16  T  ^k-1^ 

xli\  *2J/  X1X2T1  J  m/2  *  - 
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+  (P  ”  )aPX  X  T  -k  (1  s  *)6  T  ^k- 

xli\  x2j/  Wl  nU,s2  'V  2 


This  updating  breaks  down  into  particular  cases  as  follows 


1.  If  at  t,  ,  ,  x, .  =  0 
k- 1  li 


PX1X2T1T2(tk)  "  PX1X2T1T2(tk-l) 


2.  If  at  t,  x..  =  1  and  x0.  =  1 
k-l  li  2j 


P  ft  ^  «  P  (t  ) 

X1X2T1+j6m  T2v  V  X1X2T1T2'‘  k-l' 


3.  If  at  t.  x,.  =  1  and  x».  ■  0 
k-l*  li  2j 


PX1X2T1+j6m  T2(tk}  (1  "  a)PX1X2T1T2(tk-l) 


\,C  T  ^1,)  “  V  "P  T 


Vl’AU'^VV  k  WlY  k'1' 


A  combination  of  these  terms  gives 


?XlX2TlT2tk>  "  ^  "  S,lu)PxiX2TlT2<tk-1> 


+  Sx11(1  -  “(l  -  Sx2J))PX1X2I1-jJmil2(Ck-l> 

"  Si"!1  '  6l‘21)PxiX2Trkn(J’s2*)«m.I2<tk- 


e  v.*v 


y 


••  ^  a^r*  -V'**!**!?  j*"ft** 
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In  a  similar  manner,  if  at  t^  member  j  of  side  II  shoots  (according 
to  the  geometry)  at  member  i  of  side  I,  we  have  the  following  updating 
formula: 


’XlX2TlT2<tk>  '  I"  '  %j)PxiX2IlI2(tk-1> 

+  6x2.(1  -  6(1  -  Sx11))PX1X2T1I2-i6n/tk-l) 


+ 


Vli 


)PxiX2TiT2-k»(J-sl*)Jo,<tk-l) 


(3.2) 


These  launch  formulae  allow  for  the  possibility  of  continued  missile 
launchings  by  one  side  even  after  all  members  of  the  other  side  are  killed. 
This  does  not  affect  survival  results;  however,  it  is  not  realistic  in 
terms  of  missile  usage.  More  complex  updating  formulae  will  be  introduced 
for  missile  usage  in  Section  4. 

Suppose  that,  from  the  geometry  run,  t^  is  a  time  point  at  which  the 
missile  of  element  i  of  side  I  is  to  impact  some  member  of  side  II.  The 
original  geometry  would  have  specified  the  member  of  side  II  receiving 
the  impact.  Previous  impacts  on  this  particular  member  of  side  II  and  the 
assumption  of  optimistic  firings  lead  to  the  possibility  that  any  member 
of  side  II  may  be  the  target.  With  this  in  mind,  we  consider  the  fol¬ 
lowing  updating: 


P  (t  )  =  P  (t  ) 

XXX  X  k  X  X  X  T  '  k,—  ^  * 


n 


+  Yj  6X  (1  -  Pli(tk))PV  V  T  T  (tu_l) 


j=l 


2  j 


Wl^m.V 


t1+j«  i/w 

j=1\  2j /  1  2  n,  1  J  m,  2 


‘j  “i 


18 


*9*. 


NWC  TP  5815 


where  t^i»  the  ith  component  of  Tj_,  is  0.  The  above  updating  formula 
may  be  condensed  to  yield 


Px  X  T  T  -  X  T  T  ^  tk-l^ 


+  '  %JPll<tk,)PxiX2Tl+JVT2(tk-l) 


It 

+  Sf1  "  6X?  )Pli(tk)PX  X  +6  T  +j6 

j=1\  2j;  12  1J  m± 


T2(tk-1)  (3,3) 


where  tjj  =  0. 


In  a  similar  manner,  if  the  missile  fired  by  member  j  of  side  II  is 
to  impact  side  I  at  t^,  then 


Px  x  T  T  ^k^  Px  X  T  T  ^  ^'k-l^ 


+  '  6j‘liP2J(tk>)PxiX2IlI2+lln1<tk-l) 

1=1  j 


+  Sf1  "  6x  )P2j(t  )PX  4  X  T  T  4  (tk-l} 
i_1\  li /  J  k'  1  mi  2  1  2  n^ 


(3.4) 


where  -  0. 
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The  updating  for  impacts  is  not  nearly  so  straightforward  as  in 
launchings.  Consider,  for  example,  a  2  on  2  updating  where  member  1  of 
side  I  has  launched  at  side  II  with  an  impact  updating  at  t^.  Letting 
Tl  =  T2  =  O2  at  t^  and  assuming  we  want  q^q  q  (t^) ,  we  have 

PX1(l,0)0202(tk)  =  PX1(l,0)(0,0)(0,0)(tk-l) 


+  (1  Pll(tk))PX1(l,0)(l,0)02(tk-l) 


+  PX1(l,0)(2,0)02(tk-l) 


+  Pll(tk)PX1(l,l)(2,0)02(tk-l) 


A  number  of  events  may  contribute  to  the  calculation  under  consideration 
and  all  events  must  be  accounted  for  properly. 

If  at  tN  the  last  possible  impact  occurs,  we  obtain  x  0  0 

which  we  redefine  as  Pv  „  ,  1  2  m  n 

xlx2 

The  joint  probability  of  survival  of  member  i-^,  i2»  •••,  ir  of  side  I 
and  ji,  J2>  •••»  js  of  side  II  is  given  by 


K' x 


li2’  ***  Xlir’  X2j1*  ■**  X2j 


0- 


1  1 


111 


£  £  •••  £  £  £ 


xir°  x12=0  xlm'°  x2r°  x22'°  x2rT° 


{  •  •  •  X  X-.  •  •  •  X-  P 

U1  ur  21  2js  X1X2 


In  particular,  P(xn)  Is  the  probability  of  survival  of  the  ith  member  of 
side  I,  and  P(x2j)  is  the  probability  of  survival  of  the  jth  member  of 
side  II. 
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Section  4.  WEAPONS  USAGE  FOR  m  ON  n  SHOOT-LOOK-SHOOT  DUELS 


There  is  a  straightforward  way  of  expanding  the  state  space  to  keep 
track  of  missile  usage.  One  merely  lets  the  state  space  be  determined 
by  the  vector  (X^,X2 ,T^ ,T2 ,M^ ,M2)  where  Xi,X2,T]_,  and  T2  are  defined  in 
Section  3  and  is  an  m- vector  where  Mu,  the  ith  component  of  M^,  is 
the  number  of  missiles  combatant  i  on  side  I  has  used.  Similarly,  M2 
is  an  n-vector  associated  with  missile  usage  of  side  II. 

If  Li  is  the  missile  loadout  for  elements  of  side  I  and  if  L>2  is  the 
missile  loadout  for  elements  of  side  II,  then  the  state  space  has 
2m+n(n+l)ra(m+l)n(Li+l)m(L2+l)n  states.  If,  for  example,  m=n=2  and 
Li=L2=4  there  are  810,000  states  in  the  state  space. 

The  numerical  example  in  Appendix  B  suggests  that  the  situation  may 
not  be  as  bad  as  one  might  expect.  If  one  only  introduces  states  as 
required,  the  actual  number  of  states  considered  may  not  be  excessive 
at  all.  If  it  should  turn  out  that  state  space  size  starts  to  be  a 
problem,  then  Monte  Carlo  procedures  should  be  considered  for  looking 
at  this  problem  for  moderate  m,  n,  and  weapons  loadout.  The  Monte  Carlo 
approach  to  this  problem  will  be  considered  in  a  later  report. 


4.1  UPDATING  FORMULAE 

For  weapons  usage  formulae  in  this  and  the  next  section  we  have 
occasion  to  use  the  standard  Kronecker  delta. 


Also , 


mli 


=  1 

=  0 


if  b. ,  =  L. 
li  1 


iC  mu  +  Lx 


6=1 

m2j 


if  m2j  =  L2 


if  m2j  4  L2 


The  updating  formulae  are  straightforward  extensions  of  Equations  3.1, 
3.2,  3.3,  and  3.4  in  Section  3.  To  obtain  ^Xi X2T1T2M1M2 ^  ^or  anX  *■» 
we  assume 


LL  0  0  0  0 
m  n  m  n  m  n 


(t0)  =  1 

L 

where  tQ  <  t^.  Nov;  assume  PX1X2T1T2M1M2 (tk-l^  is  known. 
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If,  at  tfc,  member  i  of  side  I  shoots,  according  to  the  geometry,  at 
member  j  of  side  II,  the  updating  formula  (Equation  3.1)  is  modified 
as  follows: 

PX1X2T1T2M1M2tk)  (X  " 

+  .T2HA.M2<tk-1> 

J  11 

+  6x1,a(1  “  <^x»,)PX1X„T1  -k  (j,s*)5  T„M, -6,  .  *.5  M  (tk-l 

li  \  21’  12  1  n  J  2  m.  2  1  k  (j,s0)  m.  2 

j  j_  n  2  l 


(4.1) 


One  may  modify  Equation  3.2  in  a  similar  way  to  obtain  a  corresponding 
formula  when  member  j  of  side  II  shoots  at  member  i  of  side  I. 

If,  at  tfc,  a  missile  of  element  i  of  side  I  is  to  impact  some 
member  of  side  II,  we  modify  Equation  3.3  to  obtain 


P  It  )  =  p  (t  ) 

X1X2T1T2M1M2  k  X1X2T1T2M1M2  k_1 


+  '  {-2jPu(t,',)PxiX2Il^VI2MlM2(tk-l) 


2  26m  f1  "  6x2.)Pli(tk)PX  X  +5  T  +j6  T  M  M  -(L2~r)6  (tk-l} 
„  ,  ,  2i\  2 i/  1  2  n.  1  m.  2  1  L  L  m. 


r=0  j=l 


(4.2) 


where  t^  =  0.  A  similar  modification  of  Equation  3.4  is  made  when  a 
missile  of  element  j  of  side  II  is  to  impact  some  member  of  side  I. 


4.2  JOINT  PROBABILITY  DISTRIBUTIONS  OF  SURVIVORS  AND 
MISSILES  USED  IN  AN  m/a  ON  n/0  DUEL 

Finally  at  t^,  =  T2  =  0  and  one  obtains  the  j°int 

probability  distribution  of  survivors  and  missile  usage.  Various 
statistical  quantities  of  interest  are 
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Section  5.  SHOOT-SWITCH  FIRING  DOCTRINES 


In  this  section  we  consider  firing  doctrines  that  allow  a  combatant 
to  switch  and  fire  at  members  of  the  opposing  side  while  his  missiles 
are  in  flight.  The  only  restriction  is  that  combatant  i  of  one  side 
cannot  fire  a  missile  at  combatant  j  of  the  other  side  if  he  has  a 
missile  in  flight  at  j .  Providing  £  n,  any  member  of  side  I  can  have 
as  many  as  n  missiles  in  flight  at  a  time.  Providing  L2  £  m,  any  member 
of  side  II  can  have  as  many  as  m  missiles  in  flight  at  a  time. 

The  vectors  Xj_,  X2,  M]_,  and  M2  are  defined  as  in  previous  sections. 

We  replace  vectors  T^  and  T2  by  the  m  by  n  matrix,  Tj_,  and  the  n  by  m 
matrix,  T2.  The  entry  in  the  ath  row  and  6th  column  of  matrix  T^, 
t^ag,  is  1  if  combatant  a  on  side  k  has  a  missile  in  flight  at  combatant  6 
on  the  other  side.  This  entry  is  0  otherwise. 

With  this  notation  the  size  of  the  state  space  is  2m+n+2mn(L]+l)m 
(L2+I)11.  This  is  larger  thari  the  state  space  of  Section  4  by  a  factor  of 

**“/(.«,  W  -  (£if  (£)"  *  1 

with  exponential  growth  in  m  and  n. 

For  m=n=l,  the  ratio  is  1;  for  m=n~2,  the  ratio  is  3.16;  and  for 
m=n=3,  the  ratio  is  64.  Thus,  the  potential  size  of  the  state  space 
can  increase  dramatically  for  shoot-switch  firing  doctrines.  Appendixes  A 
and  B  indicate  that  the  number  of  states  with  non-zero  probabilities  may 
be  relatively  small  for  shoot-look-shoot  firing  doctrines.  In  fact, 
those  numerical  examples  indicate  that  the  total  number  of  states  that 
occur  with  non-zero  probabilities  in  the  entire  probability  updating  is 
relatively  small. 

5.1  BOUNDS  ON  THE  NUMBER  OF  NON-ZERO  PROBABILITY  STATES* 

We  assume  that  a  duel  consists  of  at  most  2NS  events,  the  NS  launches 
and  NS  impacts  for  all  participants.  It  is  also  assumed  that,  at  each 
launch  event,  each  state  with  a  non-zero  probability  gives  rise  to  exactly 
one  state  with  a  non-zero  probability.  Generally  speaking,  and  are 
updated  one  state  for  one  depending  on  X1X2T1T2M1M2  at  the  launch  time. 

*  This  section  is  based  on  an  analysis  of  the  state  space  size 
problem  by  Dr.  William  Alltop  of  the  Mathematical  Services  Branch,  NWC. 
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Since  the  growth  of  non-zero  probabilities  can  only  occur  at  impact 
events,  the  nature  of  this  growth  can  be  determined  by  considering  a 
single  component  i  of  Xi.  Suppose  i  of  side  I  is  impacted  Si  times  and 
all  other  components  of  Xi  and  X2  are  1.  The  vector  Xi  with  ith  com¬ 
ponent  equal  to  1  (not  killed)  is  a  possible  state  on  each  impact.  This 
same  vector  with  ith  component  equal  to  0  (killed)  may  occur  for  the  first 
time  on  any  of  the  Si  impacts.  In  general,  the  impact  number  associated 
with  the  first  time  the  zero  occurs  distinguishes  different  non-zero 
probability  states.  Consequently  the  vector  (X]_,X2)  with  a  0  in  the  ith 
component  of  Xi  and  l's  elsewhere  is  associated  with  Si  non-zero  proba¬ 
bility  states  each  having  different  entries  in  T1T2M1M2.  That  is,  for 
(Xl,X2)  as  described  there  are  Si  different  states  (X1X2T1T2M1M2)  identi¬ 
fiable  in  the  total  state  space. 


fc 

More  generally,  if  Sj  is  the  number  of  times  participant  j  on  side  II 
is  impacted  and  it  can  be  assumed  that  the  appropriate  tiij  and  t2ji  are 
all  1  for  all  impacts,  then  the  number  of  states  with  non-zero  probabili¬ 
ties  associated  with  each  (Xi,X2)  is 


m 


n  s 

i=l 


(1-6 

i 


xii' 


h  s*(1"6x2^ 

j=i  J 


A  simple  inductive  argument  shows  that 


1 


m 


n 

i=l 


Si 


n 


j=i 


*(l-<SX2j) 

sj 


m  n 

/7  d  +  si)  fj  U  +  S*) 
i=l  j=l 


When  tiij  or  t2ji  is  0,  the  possibility  for  growth  does  not  exist. 
Hence,  if  |Na|  is  the  number  of  non-zero  probability  states  at  the  final 
updating. 


m  n 

|nsI  *  77  u  +  s±)  77  (i  +  s*) 

i=i  j=i 


One  may  show  that  |NS|  is  maximized  when  impacts  are  spread  as  evenly  as 
possible;  hence, 


+ 


NS 

mfn 
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If  L]_  =  L2  =  4,  m  =  n  =  2 ,  NS  =  mL^  +  nL2  =  16, 
I  Ng |  <  625 


considerably  less  than  2.5  million.  We  now  turn  to  updating  formulae  for 
a  number  of  cases. 


5.2  SHOOT-SWITCH  WITH  a  =  0,  0  =  0 

With  a  =  3  =  0,  every  missile  scheduled  for  firing  is  fired,  if  the 
attacker  is  alive.  This  occurs  regardless  of  whether  the  target  is  alive 
or  dead.  In  this  case  the  duel  consists  of  at  most  mLi  +  nL.2  launches 
and  impacts.  The  only  way  a  participant  can  end  up  with  missiles  is  if 
he  is  not  scheduled  to  fire  all  missiles  and  he  survives  the  duel. 

The  bound  in  this  case  is 


The  comments  about  the  t^j  or  1 2 j i  at  impact  apply  here.  Also,  the  fact 
that  impact  numbers  will  not  be  exactly  mLi/n  or  nL.2/m  implies  that  this 
bound  could  be  very  conservative. 

ag 

For  this  section  we  define  6rs  to  be  an  r  by  s  matrix  of  0's  with 

the  exception  that  the  element  in  the  ctth  row  and  3th  column  is  a  1. 
Assuming  that  the  geometry  is  such  that  no  combatant  is  allowed  to  fire 
more  missiles  than  his  missile  load,  we  arrive  at  updating  formulae 
similar  to  those  in  Section  4  when  a  =  3  *  0. 

If,  at  tfe,  i  of  side  I  shoots  at  j  of  side  II,  then 


X1X2T1T2M1M2 


<ck> 


i1  -  \t)\ 


X2TlT2MlM2(tk_l) 


+  6x1,PX1X9T1-6ijT_M14  M  (tk-l) 
li  12  1  mn  2  1  m^ 


(5.1) 
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If,  at  t^,  i  of  side  I  impacts  j  of  side  II, 


P  Ct  )  =  P  (t  ) 

X1X2T1T2M1M2  k  X1X2T1T2M1M2  k-1 


+  f1  "  6x„.Pli(tk))PX  X  T  +6iJT  M  M  (tk-l) 
\  2j  /  121  nn  2  1  2 


+  i1  ~  6x  )6m  Pli(tk)  2  PX  X  +6  T  +6ijT  M  M  -(L2-k)6  (tk-l} 
\  2j/  2j  .  n  12  j  1  mn  2  1  2  2  j 


(5.2) 


where  t, , .=0. 

lij 


Similar  expressions  exist  when  j  of  side  II  launches  at  i  of  side  I 
and  when  i  of  side  I  incurs  an  impact  of  j  of  side  II. 


5.3  MODIFIED  SHOOT-SWITCH,  a  =  0,  0  =  0 

By  a  modified  shoot-switch  firing  doctrine  with  a  =  0,  3  =  0  we 
mean  one  that  allows  a  combatant  to  save  his  missile  if  his  opponent  in 
the  firing  schedule  is  dead.  This  does  not  compensate  for  the  ghost 
problem,  but  it  does  allow  one  to  conserve  missiles.  As  one  might 
expect,  the  updating  is  more  complex  and  the  duel  does  not  necessarily 
end  with  the  (mL^  +  nL^th  impact.  In  fact,  with  non-zero  probab:  1  *  ty , 
each  impact  implied  by  the  geometry  will  give  rise  to  some  new  non-zero 
(perhaps  very  small)  probability  updating. 

Letting  Sj  be  the  total  number  of  possible  impacts  on  side  I  and 
Sii  be  the  total  number  of  possible  impacts  on  side  II,  the  bound  on  the 
number  of  non-zero  probability  states  is 


Na|  £  (1  +  SI/m)m  (1  +  Sn/n)n 


If,  for  example,  Sj  =  Sn  =  16  and  m  =  n  =  2,  then  |  Ns  j  <  6561. 
However,  Sj  =  Sji  =  40  and  m  *  n  =  2  yield  |NS|  <  194,481.  As  the  duel 
progresses,  more  and  more  tijj  and  t2ji  will  be  0,  reducing  this  bound 
considerably.  One  procedure  is  to  set  aside  a  fixed  number  of  storage 
locations  (say  10,000)  and  stop  the  updating  when  those  locations  are 
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fully  utilized  or  when  the  geometry  indicates  the  duel  has  ended.  In 
most  applications  the  joint  probability  distribution  of  (Xi,X2,Mj_,M2) 
at  such  a  stopping  point  will  differ  only  slightly  from  that  same 
distribution  at  the  conclusion  of  the  duel. 

The  updating  formulae  for  i  of  I  against  j  of  II  of  the  preceding 
section  must  be  modified  to  reflect  these  considerations.  If,  at  t^, 
i  of  I  launches  at  j  of  II,  then 


PX1X2T1T2M1M2(tk)  "  f1 


'  '  64)) PWlT2MlM2<tk-l) 


+  <Sx  6X  f1_<SmMPXXT  t  M  (tk  1>  <5.3) 

li  2j\  mli'  X1X2T1  6mnT2Ml  m.M2  k_1 


If,  at  tfc,  i  of  I  impacts  j  of  II,  then  the  updating  formula  is  identical 
to  Equation  5.2. 


5.4  SHOOT-SWITCH,  a  >  0,  0  >  0 

For  the  general  case  when  a  £  0  and  8  ;>  0,  the  problem  gets  out  of 
hand.  The  reader  should  note  that,  if  tfc  is  a  time  when  a  missile  of 

1  of  I  is  to  impact  j  of  II,  X2j  =  0,  tuj  =  0,  and  tiijk  =  1  for  some 

Jk  1  ji  k  =  1,2,* ••,  there  is  no  way  of  knowing  which  one,  if  any,  of 
these  Jk  is  the  consequence  of  an  optimistic  firing  when  i  was  to  fire 
at  j.  This  presents  a  real  problem  on  how  to  update  on  impact.  It  can 
apparently  only  be  handled  by  expanding  the  state  space  to  carry  this 
type  of  information. 

One  possibility  is  to  develop  special  updating  algorithms  for 

2  on  1,  1  on  2,  and  2  on  2  duels.  Another  possibility  is  to  resort  to 
Monte  Carlo  sampling  procedures.  More  consideration  will  be  given  to 
this  problem  at  a  later  time  if  the  concepts  of  optimistic  firings  and 
shoot-switch  policies  both  have  merit  within  the  framework  of  analyzing 
aerial  combat  data. 
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Appendix  A 

A  NUMERICAL  EXAMPLE 


To  illustrate  the  use  of  the  formulae  in  Section  3,  we  shall  discuss 
in  some  detail  the  calculations  associated  with  a  sample  2  on  2  duel 
where  each  combatant  fires,  at  most,  two  missiles.  We  shall  use  the 
following  abbreviations: 

L  A  launch  event 

H  An  impact  event 


i  of  I 


Member  i  of  side  I 


j  of  II  Member  j  of  side  II 
— Launches  at 

U 

- ►  Missile  impacts 

Thus,  i  of  I  — Si-*.  j  of  II  is  read  as  member  i  of  side  I  missile  impacts 
member  j  of  side  II. 

A  typical  time  history  from  imaginary  aerial  combat  data  is  as 
follows : 


Event 


Who  against  whom 


Kill  probability  (Pfc) 


hi  =  20 


2  of  I  — 
2  of  I  — 
1  of  II  - 
1  of  II  — 

1  of  I  — 

2  of  II  — 
1  of  I  — 
1  of  II  — 

1  of  II  — 

2  of  II  — 
1  of  I  — 


1  of  II 
1  of  II 
1  of  I 

1  of  I 

2  of  II 

1  of  I 

2  of  II 
2  of  I 
2  of  I 

1  of  I 

2  of  II 
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Time 


Event 


Who  against  whom 


Kill  probability  (Pk) 


t 

t 

t 

t 

t 


12 

13 

14 

15 

16 


24 

24 

27 

28 
28 


H  1  of  I 

L  2  of  II 

L  2  of  I 

H  2  of  II 

H  2  of  I 


_L_ 

_L_ 

H 


2  of  II 
2  of  I 
2  of  II 
2  of  I 
2  of  II 


0.2 

0.7 

0.6 

0.7 

0.6 


To  illustrate  the  updating  formulae,  we  redefine  each  of  the  states 
with  which  a  nonzero  probability  is  associated  sometime  during  the  aerial 
combat.  In  fact,  states  are  only  defined  as  needed.  The  following  key 
is  required  to  understand  the  probability  updating: 


State:  Integer  State;  (Xj,  X2,  Tj,  T2) 


1 

(1,1) (1 ,1) (0,0) (0,0) 

2 

(1 ,1) (1 ,1) (0,1) (0 ,0) 

3 

(1,1) (0,1) (0,0) (0,0) 

4 

(1,1) (1,1) (0,0) (1,0) 

5 

(0,1) (1 ,1) (0,0) (0 ,0) 

6 

(1,1)  (1.1) (2 ,0) (0 ,0) 

7 

(1,1) (0,1) (2,0) (0,0) 

8 

(1»1)(1»1)(2,0)(0,1) 

9 

(0,1) (1 ,1) (0,0) (0, 1) 

10 

(0,1) (1,1) (0,0) (0,2) 

11 

(1 ,1) (0 ,1) (2 ,0) (0, 1) 

12 

(1,1)  (1 ,1) (0,0) (0,1) 

13 

(1,1) (1,0) (0,0) (0,1) 

14 

(1,1) (0,1) (0 ,0) (0,1) 

15 

(1,1) (0,0) (0,0) (0,1) 

16 

(1,1) (1 ,1) (0,0)  (2,1) 

17 

(1,1) (1 ,0) (0 ,0) (2 ,1) 

18 

(0,1) (1 ,1) (0,0) (2 ,1) 

19 

(0,1)(1,1)(0,0)(?,2) 

20 

(1 ,1) (0,1) (0,0) (0,1) 

21 

(1,1) (0,0) (0,0) (0,1) 

22 

(1 ,0) (1,1) (0,0) (0,1) 

23 

(1,0) (1,0) (0,0) (0,1) 

24 

(0,01(1,1) (0,0) (0,1) 

25 

(0,0) (1,1) (0,0) (0,2) 

26 

(1 ,0)  (1 ,1) (0,0) (0,0) 

27 

(0,0) (1 ,1) (0,0) (0,0) 

28 

(1 ,1) (1,0) (0,0) (0,0) 

29 

(0,1) (1,0) (0,0) (0,0) 

30 

(1 ,0) (1,0) (0,0)  (0,0) 

32 


.  »-•.  -t 


State:  Integer 


State:  (Xi,  X2,  Ti,  To) 


(0,0)(1,0)(0,0)(0,0) 
(0,1) (0,1) (0,0) (0,0) 

(1,1) (0,0) (0,0)  (0,0) 
(0,1) (0,0) (0,0)  (0,0) 
(1,0) (1,1) (2,0) (0,0) 

(1,1)(1,0)(2,0)(0,0) 

(1,1)(1,0)(1,0)(0,0) 
(1,0)(1,0)(2,0)(0,0) 
(1,0) (1,0) (1,0) (0,0) 

(1 ,1)  (0,0)  (2 ,0)  (0 ,0) 
(1,0) (0,0) (0,0) (0,0) 

(1,1) (1,1) (0,0) (0 ,2) 
(1,0) (1,1) (0,0) (0,2) 

(1,1) (0,1) (0 ,0) (0,2) 
(0,1) (0,1) (0,0) (0,2) 

(1,1)(1»1)(0,2)(0,2) 

(1,1) (1,0) (0,2) (0,0) 

(1,1) (1,0) (0,1) (0,0) 
(0,1) (1,1) (0 ,2) (0 ,2) 
(0,1) (1,0)  (0,2) (0,0) 
(0 ,1)  (1,0)  (0 ,1)  (0 ,0) 

(1,1) (0,1) (0,2) (0,2) 

(1,1) (0 ,0) (0,2) (0,0) 
(0,1) (0,1) (0,2) (0,2) 
(0,1) (0 ,0) (0,2) (0 ,0) 

(1,1) (1,1) (0,2) (0 ,0) 
(1 ,0) (1,1) (0,2) (0 ,0) 
(0,1) (1,1)  (0,2)  (0,0) 
(0,0) (1,1)  (0,2) (0,0) 

(1,1) (0,1) (0 ,2) (0  ,0) 
(1,0) (0,1) (0,2) (0,0) 
(0,1) (0,1)  (0,2) (0,0) 
(0,0) (0,1) (0 ,2) (0 ,0) 
(1»0) (0,1) (0,0) (0,0) 
(0,0) (0,1) (0 ,0) (0,0) 
(0,0) (0,0) (0,0) (0,0) 


-  -jr-t;  * 


i  v*r 


■ .  mu'*' 
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Nonzero 

probability 

states 

Updating  formula 

State 

probability 

Time:  tg;  action: 

•  •  • 

1 

•  •  • 

1 

Time:  ti;  action:  2  of  I  — 

_L_ 

1  of  II 

I 

I 


I 


2 

P2^1^  =  Pi^to^ 

1 

Time:  t£ ;  action:  2  of  I  - 

-JU.  II 

1 

Px(t2)  =  (1  -  Pk(t2))P2(t1) 

0.6 

3 

P3(t2)  =  Pk(t2)P2(tl) 

0.4 

Time:  t 3 ;  action:  1  of  II  — 

1  of  I 

3 

P3^fc3^  =  P3^2^ 

0.4 

4 

**4^3)  = 

0.6 

Time:  t^ ;  action:  1  of  II 

-V  1 

1 

W  *  p4(t3><1  -  pk<c4» 

0.3 

3 

P3(t4)  =  P3(t3) 

0.4 

5 

p5(t4)  -  vt3)Pk<t4) 

0.3 

Time: 


action:  1  of  1 


2  of  II 


5 

p5(t5>  =  p5(t4> 

0.3 

6 

p6(t5)  -  px<t4) 

0.3 

7 

p?(t5)  =  ?3(t4) 

0.4 

Time:  tg;  action:  2  of  II  — 

1  of  I 

8 

P8^t6)  = 

0.3 

9 

P9(t6)  =  (1  -  6)P5(t5) 

0.3(1  -  8) 

10 

Pl0(t6)  =  8P5(t5) 

0.38 

11 

pu(t6)  -  P?(t5) 

0.4 

34 
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Nonzero 

State 

probability 

probability 

states 

Updating 

formula 

Time :  t -j ; 

action:  1  of  I 

H 

- *■ 

-  11 

p9(t7)  -  PgUg) 

P10(t7)  =  P10(t6) 

P12(t7)  =  (1  '  Pk(t7))P8(t6> 

pi3(t7)  =  Pk(t7)Pg(t6) 


P14(t7> 


-  Pk(t7))Pll<t6) 


pi5(t7)  =  Pk(t7)Pu(t6) 


0.3(1 

0.38 

0.12 

0.18 

0.16 

0.24 


_ 1 

-  -  -  -  i 

Time:  tg;  action:  1  of  II  — 

u - - - —— 

2  of  I 

16 

P16(t8)  =  P12(t7) 

0.12 

17 

P17(t8)  =  P13(t7) 

0.18 

18 

P18^8^  =  P9^7^ 

0.3(1  -  8) 

19 

Pl9(tg)  =  p10(t7) 

0.3B 

20 

P20^t8^  =  P14^t7^ 

0.16 

21 

P21(t8)  =  P15(t7) 

0.24 

Time:  t9;  action:  1  of  II 

-JU.I 

9 

p9(t9)  *  0-  " 

0.06(1  -  8) 

10 

P10^C9^  ^  Pfc^g^Pig^s^ 

0.068 

12 

P12^9^  =  C1  "  Pk^^1*^^ 

0.024 

13 

Pl3(t9)  (1  “  Pk^t9))P17^t8) 

0.036 

14 

P14(t9)  =  P20(t8) 

0.16 

21 

P21(t9)  =  P21(t8) 

0.24 

22 

P22(V  "  ^^9^16^8^ 

0.096 

23 

P23(t9)  =  Pk(t9)P17(t8) 

0.144 

24 

P24(t9)  =  ^*9^18^8* 

0.24(1  -  8) 

25 

p25(t9)  =  Pk^t9)P19(t8) 

0.248 

Nonzero 
probabili  ty 
states 
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Updating  formula 


t iQ ;  action:  2  of  II 


State 

probability 


| 

Pl(t10)  (1  "  Pk(t10))P12(t9J 

0.0168 

1 

P3(t10')  (1  "  Pk(t10))P14<:t:9) 

0.112 

P5(t10j  =  Pk(t10)P12(t9) 

1 

+  P9^9^  +  ^1-Pk^t10^P10^t9^ 

0.0672  -  0.0188 

P26(t]0)  =  (1_Pk(tI0^P22(t9) 

0.0672 

P27(t10)  =  Pk(t10)P22(t9^ 

+  P24^9^  +  P2  5  ^ C  9  ^ 

i 

+  Pk(t10)P10(t9) 

=  0.3(0.096)  +  0.24(1  -  3) 

+  0.243  +  0.3(0.06)8 

0.2688  +  0.0138 

P28^t10)  ^  "  Pk^t10^P13^t9^ 

0.0252 

P29(t10)  =  Pk(t10)P13(t9) 

0.0108 

?30(t10)  (1  “  Pk(t10))P23(t9) 

0.1008 

P3l/t10'>  *  Pk^t10^P23^t9^ 

0.0432 

P32(t10}  =  Pk/t10)P14(t9) 

0.048 

P33(t10J  (1  "  Pk(t10))P2I(V 

0.168 

P34(t10J  =  Pk(,t10)P21(C9^ 

0.072 

Nonzero 

probability 

states 

Updating  formula 

State 

probability 

■immiMPi 

■2  of  11  (contd.) 

35 

P35(tlP  ~  P26(t10) 

0.0672 

36 

P36('tll^  ^  ~  a^P28^t10^ 

(1  -  a) (0.0252) 

37 

P37(tll}  =  aP28(t10) 

a(0.0252) 

38 

P38^tll^  ^  "  “^30^10^ 

(1  -  a)  (0.1008) 

39 

P39(tll)  =  aP30(t10) 

a(0.1008) 

40 

P40(tll)  (1  "  “^33**10* 

(1  -  a) (0.168) 

Time:  tj^;  action:  1  of  I 

JUn 

1 

Pl(t12>  ’  (1  -  Pk<t12))P6<Cll> 

0.01344 

3 

P3(C12)  ^  “  Pk(t12))P7(tll) 

0.0896 

5 

P5(t12)  =  P5(tll) 

0.0672  -  0.018B 

26 

P26(t12)  (1  “  Pk(t12))P35(tll) 

0.05376 

27 

P27(t12)  =  P27(tll) 

0.2688  +  0.0188 

28 

p28<ti2)  =  P6(tll)Pk(t12) 

+  P36(tll) 

+  P37(t11)(l  ~  ^^12^ 

0.02856  -  0.00504a 

29 

p29(ti2)  =  P29^tll) 

0.0108 

30 

P30(t12)  =  P35(tll)Pk(t12) 

+  P38(tll) 

+  P39(tll)(1  “  Pk(t12^ 

0.11424  -  0.02016a 

31 

P31(t12)  =  P31^tll) 

0.0432 

32 

P32(t12)  =  P32(tll) 

0.048 

33 

P33(t12)  =  Pk(t12)P7(tll) 

+  P33(tl;L)  +  P40Ct1:L) 

+  pk(t12)P37(tll) 

0.1904  +  0.00504a 

34 

P34(t12)  =  P34(tll) 

0.072 

41 

P41(t12)  =  Pk(t12)P39(tll) 

0.02016a 

0.1904  +  0.00504a 

0.072 

0.02016a 


.„!.  ,  "  ''  '.  ii  i^lr^.iliulMI»l*IIWIi>l!*ffgHr 
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Nonzero 

probability 

states 

Updating  formula 

| 

State 

probability 

— 2  of  I 

10 

P10(t13)  =  P5^t12) 

0.0672  -  0.0188 

22 

p22<ti3)  =  eP26^t12) 

8(0.05376) 

25 

P25^t13^  ^  ^P27^t12^ 

(1  -  6)(0.2688  +  0.0186) 

27 

P27(t13)  =  ^P27(t12) 

6(0.2688  +  0.0186) 

28 

P28(t13)  =  P28^t12) 

0.02856  -  0.00504a 

29 

P29(t13)  =  P29(t12^ 

0.0108 

30 

P30(t13)  =  P30(t12) 

0.11424  -  0.02016a 

31 

P31(t13)  =  P31(t12) 

0.0432 

33 

P33(t13)  =  P 33**12* 

0.1904  +  0.00504a 

34 

P 34**13*  =  P34(t12) 

0.072 

41 

P41(t13)  =  P41(t12) 

0.02016a 

42 

P42(t13)  =  Pl(t12) 

0.01344 

43 

P43(t13}  =  (1  “  S)P26(t12) 

(1  -  6)0.05376 

44 

P44(t13)  =  P3(t12) 

0.0896 

45 

P45(t13)  =  P32(t12) 

0.048 

Time:  t^4 ;  action:  2  of  I  — 

-L*.  2  of  II 

22 

P22(t14)  =  P22(t13) 

0.053766 

25 

P25(t14)  =  P25**13* 

(1  -  6) (0.2688  +  0.0186) 

27 

P27(t14)  =  P27**13* 

6(0.2688  +  0.0186) 

P30(t14)  =  P 30**13* 

0.11424  -  0.02016a 

31 

P31**14*  =  P31**13* 

0.0432 

33 

P33**14*  =  aP33**13* 

a(0.1904  +  0.00504a) 

34 

P34**14*  =  aP 34**13* 

0.072a 

41 

P41**14*  =  P41**13* 

0.02016a 

43 

P43**14*  =  P43**13* 

(1  -  6) (0.05376) 

46 

P46**14*  =  P42**13* 

0.01344 

NWC  TP  5815 


Nonzero 
probability 
s  tates 


Updating  formula 


P47^14* 

P48(t14) 

P49(t14) 

P50(t14> 

P51(t14) 

P52^t14) 

P53(t14> 

P54(tl4> 

P55(t14> 


action:  2  of  I 


(1  -  a)P2g(t13) 

aP28^13^ 

P  10^13) 

(1  -  a)P2g(t13) 

aP29 ( ^13^ 

P44(t13) 

(1  -  a)P33(t13) 

P45(t13) 

(!  -  P34 ^ ci3^ 


probability 


2  of  II  (contd.) 


(1  -  a) [0.02856  -  0.00504a] 
a[0. 02856  -  0.00504a] 

0.0672  -  0.0188 
(1  -  a) (0.0108) 
a(0 .0108) 

0.0896 

(1  -  a) [0.1904  +  0.00504a] 
0.048 

(1  -  a) (0 .072) 


e:  1 3^5 ;  action:  2  of 


26 

P26(t15)  =  P43(t14) 

+  (1  -  pk(t15))p22(ti4) 

0.05376  -  0.0376328 

27 

P27(t15}  =  Pk(t15)P22(t14) 

+  P25(t14)  +  P27(t14) 

0.2688  +  0.0556328 

30 

P30(t15)  =  P30(t14) 

0.11424  -  0.02016a 

31 

P31(t15)  =  P31(t14) 

0.0432 

33 

P33(t15)  =  P33(t14) 

a(0. 1904  +  0.00504a) 

34 

P34(t15)  =  P34(t14) 

0.072a 

41 

P41(t15)  =  P41(t14) 

0.02016a 

47 

P47(t15)  =  P47(t14) 

(1  -  a) [0.02856  -  0.00504a] 

48 

^8^15*  =  P48(t14) 

a[0. 02856  -  0.00504a] 

50 

P50^t15)  =  P50(t14) 

(1  -  a) (0.0108) 

51 

p5l(tl5)  =  P51(t1A) 

a(0 .0108) 

53 

P53(t15)  =  P53(t14) 

(1  -  a) [0.1904  +  C. 00504a] 

55 

P55(t15)  =  P55(t14) 

(1  -  a) (0.072) 

Nonzero 

probability 

states 


Updating  formula 


probability 


J _ L 


Time:  ti5;  action:  2  of  II  - ►  I  (contd.) 


56 

P56(t15)  (1  ”  Pk(t15))P46(t14) 

0.004032 

57 

P57(t15)  =  Pk(t15)P46(t14) 

0.009408 

58 

P58(t15)  (1  " 

0.02016  -  0.0054B 

59 

P59(t15)  =  ^^15^49^14* 

0.04704  -  0.01268 

60 

P60(t15)  (1  “  Pk(t15))P52(t14) 

0.02688 

61 

^l^lS^  =  Pk(t15)P52(t14) 

0.06272 

62 

P62(t15>  '  (1  -  Pk<*15»P54<V 

0.0144 

63 

_ 1 

P63(t15)  =  Pk(t15)P54(t14) 

0.0336 

tl 

Time:  t^g;  action:  2  of  I - ►  II 


1 

Pl<t16>  *  (1  -  Pk(t16))P56<t15> 

0.0016128 

3 

P3(t16)  (1  '  Pk(t16))P60(t15) 

0.010752 

5 

P5(t16>  '  (1  -  Pk<t16»P58(t15) 

0.008064  -  0.002168 

26 

P26(t16>  *  (1  -  Pk<t16»P57<t15> 

+  P26(t15) 

0.0575232  -  0.0376328 

27 

P27(t16)  (1  "  Pk(t16))P59(t15) 

+  P27(t15> 

0.287616  +  0.0505928 

28 

P28^16^  =  Pk(t16)P56(t15) 

+  P47(t15) 

0.0309792  -  0.022176a 

2 

+  p48(t15xi  -  pk<t16» 

+  0.003024a 

29 

P29(t16)  =  Pk(t16)P58(t15) 

+  P50(t15) 

0.022896  -  0.003248 

+  -  Pk(t16» 

-  0.00648a 

30 

P30(t16)  =  Pk(t16)P57(t16) 

- 

+  P30(t15> 

0.1198848  -  0.02016a 

31 

P31(t16)  = 

+  P31(t15> 

0.071424  -  0.007568 

40 


Nonzero 

probability 

states 

Updating  formula 

State 

probability 

Time:  t^g;  action:  2  of  1  - 

-►  II  (contd.) 

32 

P32(t16)  (1  "  Pk(t16))P62(t15) 

0.00576 

33 

P33('t16)  =  P33(t15)  +  P53(t15) 

0.2065280  +  0.022176a 

+  Pk(t16)[P48(t15)  +  P60(t15)] 

-  0.003024a2 

34 

P34(t16)  =  P34(t15)  +  P55(t15) 

+  Pk(t16^P51(t15)  +  P62(t15)] 

0.08064  +  0.00648a 

41 

P41(t16)  =  Pk(t16)P61(t15) 

+  r4i<‘l5> 

0.037632  +  0.02016a 

64 

P64(t16)  (1  “  Pk(t16))P61(t15) 

0.025088 

65 

P65(t16)  (1  “  Pk(t16))P63(t15) 

0.01344 

66 

P66(t16)  =  Pk(t16)P63(t15) 

0.02016 

The  survivability  of  various  participants  in  the  duel  is  given  as 
follows: 

P(x u  -  1)  =  0.49  -  0.0376326 

P(x12  *  1)  =  0.367232  -  0.00546 

P(x21  =  1)  =  0.6  -  0.048816a  +  0.003024a2 

P(x22  =  1)  =  0.409856  +  0.01086 

Also  side  I's  exchange  ratio  which  depends  on  a  and  6  is  given  by 


8(a,  6) 


2  -  (1.009856  -  0.048816a  +  0.003024a^  +  0.01086) 
2  -  (0.857232  -  0.0430326) 


0.990144  +  0.048816a  -  0.003024a4  -  0.01086 
1.142768  +  0.0430326 


where 


0.8258931  =  9(0,1)  <;  9  £  9(1,0)  =  0.9065147 
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Appendix  B 

THE  NUMERICAL  EXAMPLE  CONTINUED  TO 
INCLUDE  WEAPONS  USAGE 

In  the  previous  example  we  introduced  the  states  to  handle  the 
updating  equations  when  2^(3)2(3)2  =  1,296  states  are  possible.  Assuming 
each  participant  in  the  duel  of  Appendix  A  has  two  missiles  we  have  a 
total  of  1,296  x  32  x  32  =  104,976  possible  states.  We  shall  now  give 
a  summarized  version  of  the  updating  formulae  when  weapons  usage  is 
considered  and  we  shall  see  that  the  66  states  required  to  handle  this 
problem  for  when  weapon  usage  is  not  considered  are  increased  only  slightly. 

We  associate  with  (Mi, M2)  the  integer  to  the  base  3  and  for 
(Mu ,Mi2)  (M21  ,M22)  let  the  corresponding  state  be  Mn3^  +  Mi232  +  M2i3 
+  M22>  Since  each  of  the  states  we  consider  may  be  thought  of  as  an 
ordered  pair  of  integers  where  the  first  integer  describes  a  state  in 
Appendix  A  and  the  second  integer  is  associated  with  a  particular  missile 
count  and  varies  from  0  to  80,  we  let  the  state  space  be  described  by 
(i,j)  where  1  4  i  £  66  and  0  <  j  £  80.  The  missile  usage  integers 
required  in  the  updating  and  the  corresponding  (Mn,Mi2)  (M21»^22)  are 
as  follows: 


0 

(0,0)(0,0) 

44 

(1,1) (2,2) 

9 

(0,1)(0,0) 

52 

(1,2) (2,1) 

12 

(0,1)(1,0) 

53 

(1,2) (2,2) 

15 

(0,1) (2,0) 

66 

(2,1) (1,0) 

39 

(1,1) (1,0) 

67 

(2 ,1) (1,1) 

40 

(1,1) (1,1) 

70 

(2 ,1) (2 ,1) 

41 

(1,1)(1,2) 

71 

(2,1) (2,2) 

42 

(1 ,1) (2,0) 

79 

(2,2) (2,1) 

43 

(1,1) (2 ,1) 

80 

(2,2) (2,2) 

Thus,  18  out  of  81  possible  missile  usage  descriptions  are  required. 


43 


/fp 


Time 


Old 

state 


New 

states 


New  probabilities  (if  required) 


(31,80) 

■ 

(32,70) 

(33,71) 

VP7,70  +  P37,71>  +  P40,71  =  0'1504  ‘  °'16296“ 

(33,44) 

P ( 33 ,44 )  =  °*168a 

34 

(34,71) 

41 

C13 

10 

(10,71) 

22 

(22,80) 

25 

(25,80) 

27 

28 

(28,71) 

29 

(29,71) 

30 

(30,80) 

31 

(31,80) 

33 

(33,71) 

0.1904  -  0.16296a 

(33,44) 

0.168a 

34 

(34,71) 

41 

(41,80) 

42 

(42,71) 

43 

(43,80) 

44 

(44,71) 

45 

(45,71) 

*14 

22 

(22,80) 

25 

(25,80) 

27 

(27,79) 

30 

(30,80) 

31 

(31,80) 

33 

(33,71) 

a(0.1904  -  0.16296a) 

(33,44) 

0.168a2 

34 

(34,71) 

41 

(41,80) 

43 

(43,80) 

46 

(46,80) 

47 

(47,80) 

48 

(48,80) 

49 

(49,80) 

50 

(50,80) 

51 

(51,80) 

52 

(52,80) 

46 
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Time 

Old 

state 

New 

states 

New  probabilities  (if  required) 

t14 

53 

(53,80) 

(1  -  ct)P33  ?1  =  (1  -  a)  (0.1904  -  0.16276a) 

(53,53) 

(1  -  a)P33  44  =  (1  -  a) 0 .168a 

54 

(54,80) 

55 

(55,80) 

t!5 

26 

(26,80) 

PkP22  80  +  P25  80  =  °*2688  "  0.2131688  -  0.01882 

27 

(27,80) 

(27,79) 

P27  ?9  =  0.26888  +  0.01882 

30 

(30,80) 

31 

(31,80) 

a(0.1904  -  0.16296a) 

33 

(33,71) 

(33,44) 

0.168a2 

34 

(34,71) 

41 

(41,80) 

47 

(47,80) 

48 

(48,80) 

50 

(50,80) 

51 

(51,80) 

53 

(53,80) 

(1  -  a) (0.1904  -  0.16296a) 

(53,53) 

(1  -  a) 0.168a 

55 

(55,80) 

56  * 

(56,80) 

57 

(57,80) 

58 

(58,80) 

59 

(59,80) 

60 

(60,80) 

61 

(61,80) 

62 

(62,80) 

63 

(63,80) 

*16 

1 

(1,80) 

3 

(3,80) 

5 

(5,80) 

26 

(26,80) 

27 

(27,79) 

P2?  ?9  =  0.26888  +  0.01882 

(27,80) 

P27  80  +  (1  "  Pk)P59  =  °’287616  "  0.2182088 

-  0.01882 

28 

(28,80) 

29 

(29,80) 

, 

30 

(30,80) 

31 

(31,80) 

47 
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T.P***>  1  r' 


Old 

state 

New 

states 

32 

(32,80) 

33 

(33,71) 

p 

33,71 

(33,44) 

P33,44 

(33,80) 

VP48 
+  0 

(33,53) 

P53,53 

34 

(34,71) 

P34, 71 

(34,80) 

P55  +  1 

41 

(41,80) 

64 

(64,80) 

65 

(65,80) 

66 

(66,80) 

New  probabilities  (if  required) 


=  a(0.1904  -  0.16296a) 
=  0.168a2 


=  (1  -  a)0.168a 
=  0.072a 


kv  51 


0.08064  -  0.06552a 


The  joint  distribution  of  survival  and  weapons  usage, 
is  given  below: 


P(l,l)(l,l)(2,2)(2,2)  ‘  °-0016128 


P(l,l)(0,l)(2,2)(2,2)  =  0,010752 


(0,1)(1,1)(2,2)(2,2) 


=  0.008064  -  0.002163 


(1*0) (1,1) (2,2) (2,2) 


=  0.0575232  -  0.0376328 


P(0,0)(l,l)(2,2)(2,l)  °-26888  +  0.01832 


(0,0) (1,1) (2,2) (2,2) 


=  0.287616  -  0.2182088  -  0.0188 


(1,1) (1,0) (2,2) (2,2) 


=  0.0309792  -  0.022176a  +  0.003024a* 


(0,1) (1,0) (2,2) (2,2) 


=  0.022896  -  0.003248  -  0.00648a 


(1,0) (1,0) (2,2) (2,2) 


=  0.1198848  -  0.02016a 


(0,0) (1,0) (2 , 2) (2 , 2) 


=  0.071424  -  0.007568 


P (0,1) (0,1) (2,2) (2,2)  “  0,00576 
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P(l,l) (0,0) (2,1) (2,2) 
^(1.1) (0,0) (1,1) (2,2) 
P(l,l)(0,0)(l,2)(2,2) 
P(l,l)(0,0)(2,2)(2,2) 
P(0,l)(0,0)(2,l)(2,2) 
P  (0 , 1)  (0 ,0)  (2 , 2)  (2 ,2) 
P(l,0)(0,0)(2,2)(2,2) 
p(l,0) (0,1)  (2,2) (2,2) 
P(0,0)(0,l)(2,2)(2,2) 
P(0,0)(0,0)(2,2)(2,2) 


=  0.1904a  -  0.16296a2 
=  0.168a2 

=  0.168a  -  0.168a2 
=  0.206528  -  0.336224a  +  0.159936a2 
=  0.072a 

=  0.08064  -  0.06552a 
=  0.037632  +  0.02016a 
=  0.025088 
=  0.01344 
=  0.02016 


From  this  distribution  one  obtains  the  following  numerical  results 
for  X^,X2,Mx,  and  M2: 


E(X  )  -  1(0.357488  -  0.0430328) 

+  2(0.249872) 

=  0.857232  -  0.0430328 

E(X2)  =  1(0.300224  -  0.048816a  -  0.01088  +  0.003024a2) 
+  2(0.354816  +  0.01088) 

=  1.009856  -  0.048816a  +  0.01088  +  0.003024a2 
E(X2)  =  1.356976  -  0.0430326 
Var(X1)  =  0.622129  +  0.0307456  -  0.00185262 

E(X2)  =  1.719488  -  0.048816a  +  0.03246  +  0.003024a2 


49 
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Vai(X2)  =  0.699679  +  0.049778a  +  0.0105876 

-  0.005467a2  +  0.001054a6  -  0.000009a4 

E(X1X2)  =  1(0.1736288  -  0.003243  -  0.02664a) 

+  2(0.1073184  -  0.0397926  -  0.022176a  +  0.003024a2) 

+  4(0.0016128) 

=  0.3947168  -  0.0828246  -  0.070992a  +  0.006048a2 

CovU^)  =  -0.470964  -  0.0486266  -  0.029145a 

+  0.003456a2  -  0.002101a6  +  0.00046562 
+  0.0001306a2 

Clearly,  Xi  and  X2  display  a  strong  negative  correlation,  as  one  would 
expect. 

For  missile  usage  we  have 

E(M1)  =  2(0. 168a2)  +  3(0. 4304a  -  0.33096a2) 

+  4(1  -  0.4304a  +  0.16296a2) 

=  4  -  0.4304a  -  0.00504a2  , 

a  strictly  decreasing  function  of  a. 

E(M2)  =  3(0.26886  +  0.01862)  +  4(1  -  0.26886  -  0.01862) 

=  4  -  0.26886  -  0.01862 

E(M2)  =  4(0. 168a2)  +  9(0. 4304a  -  0.33096a2) 

+  16(1  -  0.4304a  +  0.16296a2) 

-  16  -  3.0128a  +  0.30072a2 


50 
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Var(M1)  =  0.4304a  +  0.075156a2  +  0.004338a3 

-  0.000025a4 

E(M2)  =  9(0.26886  +  0.01832)  +  16(1  -  0.26886  -  0.01862) 

=  16  -  1.88166  -  0.12662 

Var(M2)  =  0.26886  -  0.05425362  -  0.00967763  -  0.00032464 

E(M1M2)  =  8(0. 168a2)  +  12(0.26886  +  0.01862  +  0.4304a  -  0.33096a2) 
+  16(1  -  (0.26886  +  0.01862  +  0.4304a  -  0.16296a2)) 

=  16  -  1.07520  -  0.07262  -  1.7216a  -  0.02016a2 

Cov(M1M2)  =  -0 . 115692a0  -  O.OO7747a02  -  O.OO1355a20 

-  O.OOOO91a202 

E(X1M1)  =  3(0 .072a)  +  4(0.357488  -  0.0430320  -  0.072a  +  0.168a2) 

+  6 (0.3584a  -  0.33096a2) 

+  8(0.249872  -  0.3584a  +  0.16296a2) 

=  3.428928  -  0.7888a  -  0.01008a2  -  0.1721280 

Cov(X1M1)  =  -0.419847a  -  0.00576a2  -  O.O18521a0  -  O.OOO217a20 

In  a  similar  manner,  one  may  show  that 

Cov(X2M2)  =  -0.2661510  -  O.O149202  -  O.O13122a0 

-  O.OOO878a02  +  O.OOO19403  +  O.OOOO54a202 

The  resulting  negative  correlations  are  not  surprising,  since  M  =  4, 
its  largest  value,  whenever  X  =  0,  its  smallest  value. 
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